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Abstract We review theoretical aspects of unitary Fermi gas (UFG), which has been 
realized in ultracold atom experiments. We first introduce the e expansion technique 
based on a systematic expansion in terms of the dimensionality of space. We apply 
this technique to compute the thermodynamic quantities, the quasiparticle spectrum, 
and the critical temperature of UFG. We then discuss consequences of the scale and 
conformal invariance of UFG. We prove a correspondence between primary opera- 
tors in nonrelativistic conformal field theories and energy eigenstates in a harmonic 
potential. We use this correspondence to compute energies of fermions at unitar- 
ity in a harmonic potential. The scale and conformal invariance together with the 
general coordinate invariance constrains the properties of UFG. We show the van- 
ishing bulk viscosities of UFG and derive the low-energy effective Lagrangian for 
the superfiuid UFG. Finally we propose other systems exhibiting the nonrelativistic 
scaling and conformal symmetries that can be in principle realized in ultracold atom 
experiments. 
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1 Introduction 

Interacting fermions appear in various subfields of physics. The Bardeen-Cooper- 
Schrieffer (BCS) mechanism shows that if the interaction is attractive, the Fermi 
surface is unstable toward the formation of Cooper pairs and the ground state of 
the system exhibits superconductivity or superfluidity. Such phenomena have been 
observed in metallic superconductors, superfluid 3 He, and high-!T c superconductors. 
Possibilities of superfluid nuclear matter, color superconductivity of quarks, and 
neutrino superfluidity are also discussed in literatures; some of these states might be 
important to the physics of neutron stars. 

In 2004, a new type of fermionic superfluid has been realized in ultracold atomic 
gases of 40 K and 6 Li in optical traps [1,2]. Unlike the previous examples, these sys- 
tems have a remarkable feature that the strength of the attraction between fermions 
can be arbitrarily tuned through magnetic field induced Feshbach resonances. The 
interatomic interaction at ultracold temperature is dominated by binary s-wave col- 
lisions, whose strength is characterized by the s-wave scattering length a. Across the 
Feshbach resonance, a -1 can, in principle, be tuned to any value from — °° to +°°. 
Therefore, the ultracold atomic gases provide an ideal ground for studying quantum 
physics of interacting fermions from weak coupling to strong coupling. 

In cold and dilute atomic gases, the interatomic potential is well approximated 
by a zero-range contact interaction: the potential range, ro ~ 60 «o for 40 K and 
ro ~ 30flo for 6 Li, is negligible compared to the de Broglie wavelength and the mean 
interparticle distance (n -1 / 3 ~ 5000-10000 oq). The properties of such a system are 
universal, i.e., independent of details of the interaction potential. By regarding the 
two different hyperfine states of fermionic atoms as spin-| and spin-^ fermions, the 
atomic gas reduces to a gas of spin-^ fermions interacting by the zero-range poten- 
tial with the tunable scattering length a. The question we would like to understand 
is the phase diagram of such a system as a function of the dimensionless parameter 
— °° < (akp)~ l < oo, where kp = (37t 2 n) 1 ^ is the Fermi momentum. 

The qualitative understanding of the phase diagram is provided by picture of a 
BCS-BEC crossover [3, 4, 5] (the horizontal axis in Fig. 1). When the attraction 
between fermions is weak (BCS limit where akp < and \akp\ <C 1), the system is 
a weakly interacting Fermi gas. Its ground state is superfluid by the BCS mecha- 
nism, where (loosely bound) Cooper pairs condense. On the other hand, when the 
attraction is strong (BEC limit where < akp <C 1), two fermions form a bound 
molecule and the system becomes a weakly interacting Bose gas of such molecules. 
Its ground state again exhibits superfluidity, but by the Bose-Einstein condensation 
(BEC) of the tightly bound molecules. These two regimes are smoothly connected 
without phase transitions, which implies that the ground state of the system is a 
superfluid for any (akp)~ l . Both BCS and BEC limits can be understood quantita- 
tively by using the standard perturbative expansion in terms of the small parameter 
\akp\ < 1. 

In contrast, a strongly interacting regime exists in the middle of the BCS-BEC 
crossover, where the scattering length is comparable to or exceeds the mean in- 
terparticle distance; \akp\ > 1. In particular, the limit of infinite scattering length 
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Fig. 1 Extended phase diagram of the BCS-BEC crossover in the plane of the inverse scatter- 
ing length (akfi)^ 1 and the spatial dimension d. There are four limits where the system becomes 
noninteracting; akg — > ±0 and d — > 4, 2. The system is strongly interacting in the shaded region. 



\akp \ — > °°, which is often called the unitarity limit, has attracted intense attention by 
experimentalists and theorists alike. Beside being experimentally realizable in ultra- 
cold atomic gases using the Feshbach resonance, this regime is an idealization of the 
dilute nuclear matter, where the neutron-neutron scattering length a nn — -18.5 fm 
is larger than the typical range of the nuclear force ro ~ 1 .4 fm. 

Theoretical treatments of the Fermi gas in the unitarity limit (unitary Fermi gas) 
suffer from the difficulty arising from the lack of a small expansion parameter: the 
standard perturbative expansion in terms of |aA:p | is obviously of no use. Mean- 
field type approximations, with or without fluctuations, are often adopted to obtain 
a qualitative understanding of the BCS-BEC crossover, but they are not necessarily 
controlled near the unitarity limit. Therefore an important and challenging problem 
for theorists is to establish a systematic approach to investigate the unitary Fermi 
gas. 

In Sect. 2 of this chapter we describe one of such approaches. It is based on an 
expansion over a parameter which depends on the dimensionality of space [6, 7, 8]. 
In this approach, one extends the problem to arbitrary spatial dimension d, keeping 
the scattering length infinite \akp\ —> °° (the vertical axis in Fig. 1). Then we can find 
two noninteracting limits on the d axis, which are d — 4 and d — 2. Accordingly, 
slightly below four or slightly above two spatial dimensions, the unitary Fermi gas 
becomes weakly interacting and thus a "perturbative expansion" is available. We 
show that the unitary Fermi gas near d = 4 is described by a weakly interacting 
gas of bosons and fermions, while near d — 2 it reduces to a weakly interacting 
Fermi gas. A small parameter for the perturbative expansion is e = 4 — d near four 
spatial dimensions or e = d — 2 near two spatial dimensions. After performing all 
calculations treating e or e as a small expansion parameter, results for the physical 
case of d = 3 are obtained by extrapolating the series expansions to e (e) = 1 , or more 
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appropriately, by matching the two series expansions. We apply this technique, the 
e expansion, to compute the thermodynamic quantities [6, 7, 9, 10] (Sect. 2.3), the 
quasiparticle spectrum [6, 7] (Sect. 2.4), and the critical temperature [8] (Sect. 2.5). 
The main advantage of this method is that all calculations can be done analytically; 
its drawback is that interpolations to d = 3 are needed to achieve numerical accuracy. 

Then in Sects. 3 and 4, we focus on consequences of another important charac- 
teristic of the unitary Fermi gas, namely, the scale and conformal invariance [11, 12, 
13]. We introduce the notion of nonrelativistic conformal field theories (NRCFTs) 
as theories describing nonrelativistic systems exhibiting the scaling and conformal 
symmetries. In Sect. 3.1, we describe a nonrelativistic analog of the conformal al- 
gebra, the so-called Schrodinger algebra [14, 15], and show in Sect. 3.2 that there is 
an operator-state correspondence [11]: a primary operator in NRCFT corresponds 
to an energy eigenstate of a few-particle system in a harmonic potential. The scal- 
ing dimension of the primary operator coincides with the energy eigenvalue of the 
corresponding state, divided by the oscillator frequency. We use the operator-state 
correspondence to compute the energies of two and three fermions at unitarity in a 
harmonic potential exactly (Sect. 3.3 with Appendix) and more fermions with the 
help of the e expansion (Sect. 3.4). 

The enlarged symmetries of the unitary Fermi gas also constrain its properties. 
By requiring the scale and conformal invariance and the general coordinate invari- 
ance of the hydrodynamic equations, we show that the unitary Fermi gas has the 
vanishing bulk viscosity in the normal phase [12] (Sect. 4.1). In the superfluid phase, 
two of the three bulk viscosities have to vanish while the third one is allowed to be 
nonzero. In Sect. 4.2, we derive the most general effective Lagrangian for the su- 
perfluid unitary Fermi gas that is consistent with the scale, conformal, and general 
coordinate invariance in the systematic momentum expansion [13]. To the leading 
and next-to-leading orders, there are three low-energy constants which can be com- 
puted using the e expansion. We can express various physical quantities through 
these constants. 

Finally in Sect. 5, we discuss other systems exhibiting the nonrelativistic scaling 
and conformal symmetries, to which a part of above results can be applied. Such 
systems include a mass-imbalanced Fermi gas with both two-body and three-body 
resonances [16] and Fermi gases in mixed dimensions [17]. These systems can be 
in principle realized in ultracold atom experiments. 



2 e expansion for the unitary Fermi gas 

In this section, we develop an analytical approach for the unitary Fermi gas based 
on a systematic expansion in terms of the dimensionality of space by using special 
features of four or two spatial dimensions for the zero-range and infinite scattering 
length interaction. 
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2.1 Why four and two spatial dimensions are special? 



2.1.1 Nussinovs' intuitive arguments 

The special role of four and two spatial dimensions for the zero-range and infi- 
nite scattering length interaction has been first recognized by Nussinov and Nussi- 
nov [18]. At infinite scattering length, which corresponds to a resonance at zero 
energy, the two-body wave function at a short distance r — > behaves like 

where r is the separation between two fermions with opposite spins. The first singu- 
lar term ~ 1 /r cl ~ 2 is the spherically symmetric solution to the Laplace equation in 
d spatial dimensions. Accordingly the normalization integral of the wave function 
has the form 

jdrR{rfo,^ dr ± ? , (2 ) 

which diverges at the origin r — > in higher dimensions d > 4. Therefore, in the limit 
d —¥ 4, the two-body wave function is concentrated at the origin and the fermion pair 
should behave like a point-like composite boson. This observation led Nussinov and 
Nussinov to conclude that the unitary Fermi gas at d — > 4 becomes a noninteracting 
Bose gas. 

On the other hand, the singularity in the wave function (1) disappears in the limit 
d — > 2, which means that the interaction between the two fermions also disappears. 
This can be understood intuitively from the fact that in lower dimensions d < 2, 
any attractive potential possesses at least one bound state and thus the threshold of 
the appearance of the first bound state (infinite scattering length) corresponds to the 
vanishing potential. Therefore the unitary Fermi gas at d —> 2 should reduce to a 
noninteracting Fermi gas [18]. 

The physical case, d = 3, lies midway between these two limits d — 2 and d = 4. 
It seems natural to try to develop an expansion around these two limits and then ex- 
trapolate to d = 3. For this purpose, we need to employ a field theoretical approach. 



2.1.2 Field theoretical approach 

Spin- A fermions interacting by the zero-range potential is described by the following 
Lagrangian density (here and below ft = Icq = 1): 

In d = 3, the bare coupling co is related to the physical parameter, the scattering 
length a, by 
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Fig. 2 Two-fermion scattering in vacuum in the unitarity limit. The scattering amplitude iA near 
d — 4 is expressed by the propagation of a boson with the small effective coupling g, while it 
reduces to a contact interaction with the small effective coupling g 2 near d = 2. 



1 m f dk 1 / fc 2 \ 

£k = ^-) ■ (4) 



c Ana J (2n) d 2e k \ 2m J 

In dimensional regularization, used in this section, the second term vanishes and 
therefore the unitarity limit a — > °° corresponds to cq — > °°. 

In order to understand the specialty of d = 4 and d = 2, we first study the scat- 
tering of two fermions in vacuum (jj. = 0) in general d spatial dimensions. The two- 
body scattering amplitude iA is given by the geometric series of bubble diagrams 
depicted in Fig. 2. In the unitarity limit, we obtain 

iA(p ,p) = -i T79 r i (5) 

r(i-i)(- PQ + %-io+) d/2 - 1 

which vanishes whenof — >4 and d— >2 because of the poles inF (l —3). This means 
that those dimensions correspond to the noninteracting limits and is consistent with 
the Nussinovs' arguments in Sect. 2.1.1. 

Furthermore, by expanding iA in terms of e = 4 — <f <C l,we can obtain further 
insight. The scattering amplitude to the leading order in e becomes 

iA( P0 ,p) = 2 E +0 e 2 

= fe) 2 /D(po,p) + 0(e 2 ), (6) 

where we have defined g 2 = ~£ and D(po,p) = (po — + '0 + ) The latter is 
the propagator of a particle of mass 2m. This particle is a boson, which is the point- 
like composite of two fermions. Equation (6) states that the two-fermion scattering 
near d = 4 can be thought of as occurring through the propagation of an intermediate 
boson, as depicted in Fig. 2. The effective coupling of the two fermions with the 
boson is g ~ e 1 / 2 , which becomes small near d = 4. This indicates the possibility 
to formulate a systematic perturbative expansion for the unitary Fermi gas around 
d = 4 as a weakly interacting fermions and bosons. 
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Similarly, by expanding iA in e = d — 2 <C 1, the scattering amplitude becomes 

2iz 

iA(po,p) = i—e + 0(e 2 ) = ig 2 + 0(e 2 ) 7 (7) 
m 

where we have defined g 2 =—e. Equation (7) shows that the two-fermion scattering 
near d — 2 reduces to that caused by a contact interaction with the effective coupling 
g 2 (Fig. 2). Because g 2 ~ I is small near d — 2, it will be possible to formulate 
another systematic perturbative expansion for the unitary Fermi gas around d = 2 as 
a weakly interacting fermions. 



2.2 Feynman rules and power counting of e 

The observations in Sect. 2.1.2 reveal how we should construct the systematic ex- 
pansions for the unitary Fermi gas around d — 4 and d — 2. Here we provide their 
formulations and power counting rules of e (e). The detailed derivations of the power 
counting rules can be found in Ref. [7]. 



2.2.1 Around four spatial dimensions 

In order to organize a systematic expansion around d = 4, we make a Hubbard- 
Stratonovich transformation and rewrite the Lagrangian density (3) as 



d 2 \ 1 

£^ E v4 ( idt + 2^ + ^) ^ ff " ^Q + vlvlQ+QWiYt ( 8a ) 

/ Q2 \ i 
= l F t /5 f + 03— + C73JU MP t + f t C7 + , f / + t, f t (7_ , f, (8b) 

V 2m / c 

where f = V^) r is a two-component Nambu-Gor'kov field, a± = \ (a\ ± ;<72), 
and di 2,3 are the Pauli matrices. Because the ground state of the system at finite 
density is a superfluid, we expand (j) around its vacuum expectation value 0o = 

($) >0as 



= 00 + g(P with g = ± '- — . (9) 

m \ 2% J 

Here we introduced the effective coupling g and chose the extra factor (?n(j)()/27z) e ^ 4 
so that <p has the same dimension as a nonrelativistic field. 1 



The choice of the extra factor is arbitrary, if it has the correct dimension, and does not affect final 
results because the difference can be absorbed by the redefinition of q>. The particular choice of g 
in Eq. (9) [or g in Eq. (13)] simplifies expressions for loop integrals in the intermediate steps. 



Unitary Fermi gas, e expansion, and nonrelativistic conformal field theories 



9 



Because the Lagrangian density (8b) does not have the kinetic term for the boson 
field <p, we add and subtract its kinetic term by hand and rewrite the Lagrangian 
density as a sum of three parts, £ = £q + £\ + £2, where 

£0 = f f (id, + 03^ + Oifc) ^ + <P t (&t + ^) 9, dOa) 
£1 = g^o + Y<p + g<pWo-Y + J txf /t c73f / + 2 J u(p t (p> (10b) 
£2 = -(p r (id, + ^J(p-2lKp^(p. (10c) 

Here we set 1 /cq = in the unitarity limit. We treat £0 as the unperturbed part 
and £\ as a small perturbation. Note that the chemical potential jj. is also treated 
as a perturbation because we will find it small, ji /0o ~ e, by solving the gap equa- 
tion. Physically, £0 + £1 is the Lagrangian density describing weakly interacting 
fermions f and bosons (p with a small coupling g ~ e 1 / 2 . £2 plays a role of counter 
terms that cancel 1/e singularities of loop integrals in certain types of diagrams 
(Fig. 3). 

The unperturbed part £0 generates the fermion propagator, 

g(po,p)= 2 ^ ( po i\;t), (id 

p^-E p + iO+ \ -<h Po-£ P J 



with E p = \/£p + % and the boson propagator 



D(P0,P) = ^— — ■ (12) 

Po~ i f + «0+ 

The first two terms in the perturbation part £\ generate the fermion-boson vertices 
whose coupling is g in Eq. (9). The third and fourth terms are the chemical potential 
insertions to the fermion and boson propagators. The two terms in £2 provide addi- 
tional vertices, —iU^ and — 2iji in Fig. 3(b) and (d), to the boson propagator, where 
n (po,p) =po-^-. 

The power counting rule of e is simple and summarized as follows: 

1. We regard 0o as 0(1) and hence n ~ £0 O as 0(e). 

2. We write down Feynman diagrams for the quantity of interest using the propaga- 
tors from £0 and the vertices from £ 1 . 

3. If the written Feynman diagram includes any subdiagram of the type in Fig. 3(a) 
or Fig. 3(c), we add the same Feynman diagram where the subdiagram is replaced 
by the vertex from £2 in Fig. 3(b) or Fig. 3(d). 

4. The power of e for the given Feynman diagram is 0(e N s/ 2+N f), where N g is the 
number of couplings g and Nn is the number of chemical potential insertions. 

5. The only exception is the one-loop vacuum diagram with one jj. insertion (the 
second diagram in Fig. 5), which is 0(1) instead of 0(e) due to the 1 /e singular- 
ity arising from the loop integral. 
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(b) 
-2ifi 

(d) 



0(e) 



= 0(e 2 ) 



Fig. 3 Power counting rule of e. The 1/e singularity in the boson self-energy diagram (a) or (c) can 
be canceled by combining it with the vertex from £?> (b) or (d), to achieve the simple e counting. 
Solid (dotted) lines are the fermion (boson) propagators ;'G (;'D). The fennion loop in (c) goes 
around clockwise and counterclockwise and the cross symbol represents the fl insertion. 



The same or similar power counting rule can be derived for the cases with unequal 
chemical potentials fj^ ^ ji^ [7], unequal masses ^ [19], at finite temperature 
T y^O [8], and in the vicinity of the unitarity limit cq^O [7]. 



2.2.2 Around two spatial dimensions 

The systematic expansion around d = 2 in terms of e can be also organized in a 
similar way. Starting with the Lagrangian density (8b), we expand (j) around its 
vacuum expectation value 0o = (0) > as 

= o + g<p with g=[ (-£) . (13) 



\ m J \ 2k J 

Here we introduced the effective coupling g and chose the extra factor (?n)i j2%) ~ e l A 
so that <p^<p has the same dimension as the Lagrangian density. 2 

In the unitarity limit 1/co = 0, we rewrite the Lagrangian density as a sum of 
three parts, £ = Co + C\ + £2, where 

Co = f /t (i% + cj 3 ^ + ct 3j u + (Jj foj f , (14a) 

£1 = -^V + ^c+^ + ^^^-f, (14b) 
£2 = <?V (14c) 

We treat £0 as the unperturbed part and C\ as a small perturbation. Physically, £0 + 
£1 is the Lagrangian density describing weakly interacting fermions f . Indeed, if 
we did not have £2, we could integrate out the auxiliary fields <p and (p\ 



2 See the footnote after Eq. (9). 
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Fig. 4 Power counting rule of e. The 1/e singularity in the self-energy diagram of (p field (a) can 
be canceled by combining it with the vertex from £2 (b) to achieve the simple e counting. Solid 
(dotted) lines are the fermion (auxiliary field) propagators iG (iD). 



C l ^g 2 ^a + ^a^ = g 2 wlwlWiWt, (15) 

which is the contact interaction between fermions with a small coupling g 2 ~ e. £2 
plays a role of a counter term that cancels 1/e singularities of loop integrals in a 
certain type of diagrams (Fig. 4). 

The unperturbed part £0 generates the fermion propagator 

G(po,P) = — — = 2 ■ , n+ 1 n . : , .. (16) 



with E p = J (e p — /i) 2 + . The first term in the perturbation part L\ generates 

the propagator of <p field, D(po,p) = — 1, and the last two terms generate the ver- 
tices between fermions and <p field with the coupling g in Eq. (13). Co provides an 
additional vertex i in Fig. 4(b) to the <p propagator. 

The power counting rule of I is simple and summarized as follows: 

1. We regard as 0(1). 

2. We write down Feynman diagrams for the quantity of interest using the propaga- 
tor from £0 and the vertices from C\ . 

3. If the written Feynman diagram includes any subdiagram of the type in Fig. 4(a), 
we add the same Feynman diagram where the subdiagram is replaced by the 
vertex from Co in Fig. 4(b). 

4. The power of e for the given Feynman diagram is 0(e Ns l 2 ), where Ng is the 
number of couplings g. 



2.3 Zero temperature thermodynamics 

We now apply the developed e expansion to compute various physical quantities of 
the unitary Fermi gas. Because of the absence of scales in the zero-range and infi- 
nite scattering interaction, the density n is the only scale of the unitary Fermi gas at 
zero temperature. Therefore all physical quantities are determined by simple dimen- 
sional analysis up to dimensionless constants of proportionality. Such dimensionless 
parameters are universal depending only on the dimensionality of space. 
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Fig. 5 Vacuum diagrams contributing to the effective potential near d = 4 to leading and next-to- 
leading orders in e. The boson one-loop diagram vanishes at zero temperature. 



A representative example of the universal parameters is the ground state energy 
of the unitary Fermi gas normalized by that of a noninteracting Fermi gas with the 
same density: 

£ — ^unitary (yj^ 
^free 

sometimes called the Bertsch parameter, measures how much energy is gained 
due to the attractive interaction in the unitarity limit. In terms of the pressure 
P, the energy density £, the chemical potential fi, and the sound velocity v s of the 
unitary Fermi gas are given by 




(18) 

a -t- z tp« a -f- l tp vp v a 

where £p = hp /2m and Vp = kp/m are the Fermi energy and velocity with kp = 

\2 d ~ 1 n d l 2 r ( i + l) «] being the Fermi momentum in d spatial dimensions. ^ is 
thus a fundamental quantity characterizing the zero temperature thermodynamics of 
the unitary Fermi gas. 



2.3.1 Next-to-leading orders 



4 can be determined systematically in the e expansion by computing the effective 
potential V e ff (0o)> whose minimum with respect to the order parameter <po provides 
the pressure P= min V e ff(0o)- To leading and next-to-leading orders in e, the effec- 
tive potential receives contributions from three vacuum diagrams depicted in Fig. 5. 
The third diagram, fermion loop with one boson exchange, results from the sum- 
mation of fluctuations around the classical solution and is beyond the mean field 
approximation. Any other diagrams are suppressed near d — 4 by further powers of 
e. Performing the loop integrations, we obtain 



V e ff ( 0o ) 



7-3(y+ln2) 



3Ce 



H] l-2(y-ln2) c 



i/2 



mfo 



(19) 
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Fig. 6 Vacuum diagrams contributing to the effective potential near d = 2 to leading and next-to- 
leading orders in e. 



where C pa 0.14424 is a numerical constant. The minimization of V e ff(0o) with re- 
spect to 0o gives the gap equation, d V e ff /d(j)o = 0, which is solved by 

= — [1 + (3C-1 +ln2)e + 0(e 2 )] . (20) 

The effective potential at its minimum provides the pressure P(jJ.) = — V e ff(0o(jU)) 
as a function of jj.. From Eq. (18) and n = dP/dfi, we obtain h, expanded in terms 
of e: 

l-|3C-^(l-ln2)|e + 0(e 2 ) 

= I e 3/2 + l e 5/2j ne _ 0x)246e 5/2 + 0{e V2y (21) 

2 16 

Although our formalism is based on the smallness of e, we find that the next-to- 
leading-order correction is quite small compared to the leading term even at e = 1 . 
The naive extrapolation to the physical dimension d = 3 gives 

4^0.475 (e^l), (22) 

which is already close to the result from the Monte Carlo simulation | pa 0.40( 1 ) [20] . 
For comparison, the mean field approximation yields | pa 0.591. 

The above result can be further improved by incorporating the expansion around 
d = 2. Because we can find by solving the gap equation that the order parameter is 
exponentially small [7] 

Ull)=^e- l ' l - x+0 ^\ (23) 

its contribution to the pressure is negligible compared to any powers of e. To leading 
and next-to-leading orders in e, the effective potential receives contributions from 
two vacuum diagrams depicted in Fig. 6. Any other diagrams are suppressed near 
d = 2 by further powers of e. From the pressure 

we obtain £ expanded in terms of e: 



/ 6ln £ \ 
^4 = ^-exp(-— ) 



14 



Yusuke Nishida and Dam Thanh Son 




Fig. 7 t, as a function of the spatial dimension d. Left panel: The upper (lower) curve is the ex- 
trapolation from the NLO expansion around d = 4 in Eq. (21) [d = 2 in Eq. (25)]. The middle three 
curves show the Pade interpolations of the two NLO expansions. The symbol aid = 3 indicates the 
result E, i=a 0.40(1) from the Monte Carlo simulation [20]. Right panel: The same as the left panel 
but including the NNLO corrections in Eqs. (29) and (30). 



<^2 = l-e + 0(6 2 ). (25) 

The value of E, in d = 3 can be extracted by interpolating the two expansions 
around d = 4 and d = 2. The simplest way to do so is to employ the Pade approxi- 
mants. We write Ej in the form 

e 3 / 2 / elne \ 

^=- e H^) F(e) ' <26) 

where F(e) is an unknown function and the trivial nonanalytic dependence on e 
was factored out. 3 We approximate F(e) by ratios of two polynomials (Pade ap- 
proximants) and determine their coefficients so that E, in Eq. (26) has the correct 
next-to-leading-order (NLO) expansions both around d = 4 [Eq. (21)] and d = 2 
[Eq. (25)]. The left panel of Fig. 7 shows the behavior of E, as a function of d. 
The middle three curves are the Pade interpolations of the two NLO expansions. In 
d = 3, these interpolations give 

| w 0.391, 0.366, 0.373, (27) 

which span a small interval E, rj 0.377 ± 0.014. 

We can also employ different interpolation schemes, for example, by applying 
the Borel transformation to F(e) = | fo'dt e~'l e G{t) and approximating the Borel 
transform G(t) by the Pade approximants [7]. These Borel-Pade interpolations of 
the two NLO expansions give 



3 It has been shown that F(e) has a nonanalytic term — |e 3 lne to the next-to-next-to-next-to- 
leading order in e [9]. Because we are working up to 0(e 2 ), we neglect such a nonanalytic contri- 
bution. 
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% ps 0.391, 0.364, 0.378 (28) 

in d = 3, which again span a small interval £, ps 0.378 ±0.013. Comparing the results 
in Eqs. (27) and (28), we find that the interpolated values do not depend much 
on the choice of the Pade approximants and also on the employment of the Borel 
transformation. 



2.3.2 Next-to-next-to-leading orders 

The systematic calculation of | has been carried out up to the next-to-next-to- 
leading-order (NNLO) corrections in terms of e [9]: 

^ 4 = 4- exp (^fe) [1- 0.04916e- 0.95961 e 2 + 0(e 3 )] (29) 

and in terms of e [10, 19]: 

|| d _, 2 = 1 -e + 0.80685 f + 0(e 3 ). (30) 

The middle four curves in the right panel of Fig. 7 are the Pade interpolations of the 
two NNLO expansions. In d = 3, these interpolations give 

4 ps 0.340, 0.372, 0.370, 0.357, (31) 

which span an interval t, ps 0.360 ± 0.020. 4 In spite of the large NNLO corrections 
both near d — 4 and d = 2, the interpolated values are roughly consistent with the 
previous interpolations of the NLO expansions (compare the two panels in Fig. 7). 
This indicates that the interpolated results are stable to inclusion of higher-order cor- 
rections and thus the e expansion has a certain predictive power even in the absence 
of the knowledge on the higher-order corrections. 

Finally we note that the limit % \d^,A ^ is consistent with the Nussinovs' picture 
of the unitary Fermi gas as a noninteracting Bose gas and £, \d-^i — >• 1 is consistent 
with the picture as a noninteracting Fermi gas in Sect. 2.1.1. It would be interesting 
to consider how £, should be continued down to d —¥ 1 . The unitary Fermi gas in 
d = 3 is analytically continued to spin- 1 fermions with a hard-core repulsion in d = 
1, which is equivalent for the thermodynamic quantities to free identical fermions 
with the same total density [21]. Therefore it is easy to find 

SL*-i->4. (32) 

The incorporation of this constraint on the Pade interpolations of the two NNLO 
expansions yields £ ps 0.365 ± 0.010 in d = 3. 



4 If we excluded the interpolation by the simple polynomial (| ps 0.340) as was done in Ref. [9], 
we would obtain 0.367 ± 0.010, which is consistent with the Borel-Pade interpolations without the 
NNLO correction near d = 2 [9]. 



16 



Yusuke Nishida and Dam Thanh Son 



a>03 = > X 



-iE = ^ » + ->~^ » ^->— + 0(e 2 ) 

Fig. 8 Fermion self-energy diagrams near d = 4 to the order 0(e). 



2.4 Quasiparticle spectrum 

The e expansion is also useful to compute other physical quantities of the unitary 
Fermi gas. Here we determine the spectrum of fermion quasiparticles in a superfluid. 
To leading order in e, the dispersion relation o (p) is given by the pole of the fermion 

propagator in Eq. (11); <X)(p) = E p = J £ p + <Pq- It has a minimum at \p\ = with 

the energy gap equal to A = <po = 2ji je. 

To the next-to-leading order in e, the fermion propagator receives corrections 
from three self-energy diagrams depicted in Fig. 8. We can find that E ~ 0(e) is 
diagonal with the elements 

v i \ 8 2 f dk E k -e k 

and Exi(POi p) = —En(—P0i — P)- From the pole of the fermion propagator, 
det[G _1 (ft),p) + /KT3 — E(co,p)] = 0, we obtain the dispersion relation 

»(f0 =E P + + £ "'g' 2 ^ + Ofc 2 ). (34) 

where £n and £22 are evaluated at po = E p . 

The minimum of (o(p) will appear at small momentum e p ~ 0(e). By expanding 
fi)(p) with respect to e p , we find that the dispersion relation around its minimum has 
the form 

co(p)cA+ { -^-^. (35) 
Here A is the energy gap of the fermion quasiparticle given by 
A |^4 = </> [l-(81n3-121n2)e + <9(e 2 )] 



— [l + (3C-l-81n3 + 131n2)e + 0(e 2 )] . (36) 



The minimum of the dispersion relation is located at e p = £0 > with 

eo\ d ^ 4 =ll+ e -Y+0(e 2 )=2ii + 0(e 2 ), (37) 
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where the solution of the gap equation (20) was substituted to <pQ. 

We can see that the next-to-leading-order correction is reasonably small com- 
pared to the leading term at least for A//1 = 2/e (1 — 0.345 e). The naive extrapola- 
tion to the physical dimension d = 3 gives 

H.31 and —^2 fe-M), (38) 

both of which are quite close to the results from the Monte Carlo simulation 
A/jj, rj 1.2 and Eo/fi w 1.9 [22]. For comparison, the mean field approximation 
yields 1.16 and £o//i « 1. 

Frequently A and £o are normalized by the Fermi energy £ F . By multiplying 
Eqs. (36) and (37) by fife? = ^ obtained previously in Eq. (21), we find 

— ^0.604 and — -> 1 (e->l). (39) 

£ F £f 

For comparison, the Monte Carlo simulation gives A/ef « 0.50 and £o/£f ~ 
0.8 [22] and the mean field approximation yields A /ep rj 0.686 and £q/£f ~ 0.591. 



2.5 Critical temperature 

The formulation of the e expansion can be extended to finite temperature T ^ by 
using the imaginary time prescription [8]. Here we determine the critical tempera- 
ture T c for the superfluid phase transition of the unitary Fermi gas. 

The leading contribution to T c /ef in terms of e can be obtained by the following 
simple argument. In the limit d — s- 4, the unitary Fermi gas reduces to a noninter- 
acting fermions and bosons with their chemical potentials fi and 2/i, respectively 
[see Eq. (10)]. Because the boson's chemical potential vanishes at the BEC critical 
temperature, the density at T = T c is given by 




Comparing the contributions from Fermi and Bose distributions, we can see that 
only 8 of 9 fermion pairs form the composite bosons while 1 of 9 fermion pairs is 

dissociated. With the use of £ F = ^ [r(| + l) §] 21 d at d = 4, we obtain 



7c 

£ F 



(41) 
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0(1) 0(1) 0(e) 

Fig. 9 Three types of diagrams contributing to the effective potential at finite temperature near 
d = 4 to leading and next-to-leading orders in e. Each [l (2jli) insertion to the fermion (boson) 
line increases the power of e by one. Note that at T > T c where 0o = 0, the insertion to the first 
diagram does not produce the 1 /e singularity. 



In order to compute the 0(e) correction to T c /ef, we need to evaluate the effective 
potential V e ff(0o) at finite temperature. The effective potential to leading order in e 
is given by the fermion one-loop diagrams with and without one ji insertion and the 
boson one-loop diagram. Near T ~ T c , we can expand Veff(^o) m $o/T and obtain 



Veff(0)-Veff(O)' 



rin2 



J" 



m<j)Q 
lit 



16T 



m<j)Q 
~lit 



(42) 



The second order phase transition occurs when the coefficient of the quadratic term 
in 0o vanishes. Accordingly, at T = T c , the chemical potential is found to be 



^=6^^2 + 0(6 



2, 



(43) 



We then compute the 0(e) correction to the density n = —dV e ff/djj. in Eq. (40) at 
T = T c and hence </>o = 0. For this purpose, we need to evaluate the fermion and 
boson one-loop diagrams with two ji insertions and the two-loop diagrams with 
one /i insertion to the fermion and boson propagators (Fig. 9). Performing the loop 
integrations and substituting Eq. (43), we can find 



3n 2 f37r 2 ln2- 
~2~ 



18f'(2) 



-D-2(ln2) z \e 



where Dwl.92181isa numerical constant. From the definition £p = 
we obtain T c /ep up to the next-leading order in e: 



_ In 



7c 



= 0.260-0.0112e + O(e 2 ). 



[r 



(45) 



On the other hand, because the unitary Fermi gas near d — 2 reduces to a 
weakly interacting Fermi gas, the critical temperature is provided by the usual 
BCS formula T c = (e"t /it) A, where A is the energy gap of the fermion quasipar- 
ticle at zero temperature. Because A in the expansion over e is already obtained by 
00 = 2/ie- l ^- l+0 ^ in Eq. (23), we easily find 
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Fig. 10 Critical temperature T c as a function of the spatial dimension d. The upper (lower) curve 
is the extrapolation from the NLO expansion around d — 4 in Eq. (45) [d = 2 in Eq. (46)]. The 
middle three curves show the Pade and Borel-Pade interpolations of the two NLO expansions. The 
symbol at d = 3 indicates the result T c /ep = 0.152(7) from the Monte Carlo simulation [23]. 



?f!,-i/e-i+o(«) 



(46) 



where we used n/Sp = % = 1+ 0(e) [Eq. (25)]. The exponential e~ x l e is equivalent 
to the mean field result while the correction e _1 corresponds to the Gor'kov-Melik- 
Barkhudarov correction in d = 3. 

Now the value of T c /£p in d = 3 can be extracted by interpolating the two expan- 
sions around d = 4 and d = 2 just as has been done for t, in Sect. 2.3.1. We write 
T c /£p in the form 

^ = —e- l,e - l F{l), (47) 

bp 71 

where F(e) is an unknown function and the nonanalytic dependence on e was fac- 
tored out. We approximate F(e) or its Borel transform G(t) by ratios of two poly- 
nomials and determine their coefficients so that T c /ep in Eq. (47) has the correct 
NLO expansions both around d = 4 [Eq. (45)] and d = 2 [Eq. (46)]. Fig. 10 shows 
the behavior of T c /ep as a function of d. The middle three curves are the Pade and 
Borel-Pade interpolations of the two NLO expansions. In d = 3, these interpolations 
give 

i 0.173, 0.175, 0.192, (48) 



which span an interval T c /ep ps 0.180±0.012. Our interpolated values are not too 
far from the result of the Monte Carlo simulation T c /ep = 0.152(7) [23]. For com- 
parison, the mean field approximation yields T c /bp sa 0.496. 
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2.6 Phase diagram of an imbalanced Fermi gas 



The phase diagram of a density- and mass-imbalanced Fermi gas was studied using 
the e expansion in Refs. [7, 19]. 



3 Aspects as a nonrelativistic conformal field theory 

Because of the absence of scales in the zero-range and infinite scattering interaction, 
the theory describing fermions in the unitarity limit is a nonrelativistic conformal 
field theory (NRCFT). In this section, after deriving the Schrodinger algebra and 
the operator-state correspondence in NRCFTs, we compute the scaling dimensions 
of few-body composite operators exactly or with the help of the e expansion, which 
provide the energies of a few fermions at unitarity in a harmonic potential. 



3.1 Schrodinger algebra 

We start with a brief review of the Schrodinger algebra [14, 15]. For the latter ap- 
plication, we allow spin-f and 4- fermions to have different masses nu and m^. We 
define the mass density: 

p{x)= £ m c yl(x)y a {x) (49) 

<T=U 

and the momentum density: 

ji{x) = - 1 - £ [yl{x)diY<j(x) - diYl(x)Y a (x)] , (50) 

<T=U 

where i = l,...,d and the arguments of time are suppressed. The Schrodinger alge- 
bra is formed by the following set of operators: 
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Table 1 Part of the Schrodinger algebra. The values of [X, Y] are shown below. 



X\Y 


Pj 


Kj 


D 


C 


H 


Pi 





-iSijM 


-'Pi 


-iKj 





Ki 







iK, 





(Pj 


D 


iPj 


-iKj 





-ZiC 


2iH 


C 


iKj 





2iC 





iD 


H 





-iPj 


-2M 


-iD 






mass: M = jdxp(x) (51) 
momentum : P, = J dx ji(x) (52) 
angular momentum : J,j = j dx [xjjj(x) —Xjjj(x)] (53) 



Galilean boost : Ki = dxxip(x) (54) 



H= > Mz- 

2m c 



dilatation: D = Jdxx-j(x) (55) 
special conformal : C = dx—p(x) (56) 
and the Hamiltonian: 

dyt{x) -d^aix) 

j uX 

tr=t,4-" 

+ JdxJdy\i/l(x)\i/l(y)V(\x-y\)\i/ i (y)\i/ t (x). (57) 

D and C are the generators of the scale transformation (x — >• e^a:, f — >• e 2 ^f) and the 
special conformal transformation [a; — > x/(l +Xt), t — K/(l + Xt)], respectively. 
In a scale invariant system such as fermions in the unitarity limit, these operators 
form a closed algebra. 5 

Commutation relations of the above operators are summarized in Table 1 . The 
rest of the algebra is the commutators of M, which commutes with all other op- 
erators; [M, any] = 0. The commutation relations of /y with other operators are 
determined by their transformation properties under rotations: 

[7 i7 , N] = [Jtj, D] = [Jij, C] = [Jij, H] = 0, (58a) 
[Jij, Pk] = i($ikPj - S jk Pi), [J ih K k ] = i{8 ik Kj - 8 jk K,) : (58b) 
[Jij, Jki] = i($kJji + SjiJik - dnJjk - 8 jk Ju). (58c) 



5 One potential that realizes the unitarity interaction is V(r) = (k/2) 2 lim ro ^o 6(>"o — r )/C^ m ti r o)' 
where m^i = m^m^j (m^ + m^) is the reduced mass. 
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These commutation relations can be verified by direct calculations, while only the 
commutator of [D, H] — 2iH requires the scale invariance of the Hamiltonian in 
which the interaction potential has to satisfy V(e^r) = e~ 2 *"V(r). 



3.2 Operator-state correspondence 
3.2.1 Primary operators 

We then introduce local operators 0(t,x) as operators that depend on the position 
in time and space (t,x) so that 

0(t,x) = e 'H'-iPi*iO(0)e-' Ht+,PiXi . (59) 

A local operator O is said to have a scaling dimension Aq and a mass Mq if it 
satisfies 

[D,O(0)]=iA o O(0) and [M, O(0)] =M o O(0). (60) 

With the use of the commutation relations in Table 1, we find that if a given local 
operator O has the scaling dimension Aq, then the scaling dimensions of new local 
operators [P u O], [H, O], [K h O], and [C, O] are given by Aq + 1, A a +2, Aq - 1, 
and Aq — 2, respectively. Therefore, by repeatedly taking the commutators with Kj 
and C, one can keep to lower the scaling dimensions. However, this procedure has 
to terminate because the scaling dimensions of local operators are bounded from 
below as we will show in Sect. 3.2.3. The last operator pr j obtained in this way 
must have the property 

[^,Op ri (0)] = [C,O pri (0)]=0. (61) 

Such operators that at t = and x = commute with Kj and C will be called primary 
operators. Below the primary operator pr j is simply denoted by O. 

Starting with an arbitrary primary operator 0(t,x), one can build up a tower of 
local operators by taking its commutators with P; and H, in other words, by taking 
its space and time derivatives (left tower in Fig. 11). For example, the operators with 
the scaling dimension Aq + 1 in the tower are [Pi, O] = idjO. At the next level with 
the scaling dimension Aq +2, the following operators are possible; [Pi, [Pj, O]] = 
— <9,<9,0 and [H, O] = —id t O. Commuting those operators with Kj and C, we can 
get back the operators into the lower rungs of the tower. The task of finding the 
spectrum of scaling dimensions of all local operators thus reduces to finding the 
scaling dimensions of all primary operators. 

It is worthwhile to note that two-point correlation functions of the primary op- 
erators are constrained by the scale and Galilean invariance up to an overall con- 
stant [24]. For example, the two-point correlation function of the primary operator 
O with its Hermitian conjugate is given by 
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Pi 



H 



K, 



— A Cpii +2 
C A Cpii + l 

— A o P „ 



Qi 



l*ow> 



Fig. 11 Correspondence between the spectrum of scaling dimensions of local operators in NRCFT 
(left tower) and the energy spectrum of states in a harmonic potential (right tower). The bottom of 
each tower corresponds to the primary operator pr j (primary state I'ifo^)). 



/ 2 \A -d/2-\ 

This formula or its Fourier transform I — po + jm~ ] — '0 + J w ^ t> e use " 

ful to read off the scaling dimension Aq. 



3.2.2 Correspondence to states in a harmonic potential 

We now show that each primary operator corresponds to an energy eigenstate of 
the system in a harmonic potential. The Hamiltonian of the system in a harmonic 
potential is 

H m =H + co 2 C, (63) 

where (0 is the oscillator frequency. We consider a primary operator 0(t,x) that is 
composed of annihilation operators in the quantum field theory so that (t,x) acts 
nontrivially on the vacuum: T |O) 7^ 0. We construct the following state using 
put at t = and x = 0: 

\V o }=e- H /<°O\0)\0). (64) 

If the mass of is M c t > 0, then |¥fo) is a mass eigenstate with the mass eigen- 
value M f. M\Wo) = M i -|fe>)- Furthermore, with the use of the commutation 
relations in Table 1, it is straightforward to show that \Wo) is actually an energy 
eigenstate of the Hamiltonian H m with the energy eigenvalue E = AqCO: 

Hol^o) = {H + co 2 C)e- H / m O\0)\0) 

= e- H l a {orC-icoD)o\Q)\0) = (oA o \ x I'o)i (65) 

where we used [C, O f (0)] = and [D, + (O)] = iA and the fact that both C and D 
annihilate the vacuum. 
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Table 2 Commutation relations of the Hamiltonian and the raising and lowering operators in the 
oscillator space. The values of [X, Y] are shown below. 



X\Y 




Qj 


H a /w 


L 


V 


Ql 





-SijM 


-Ql 


-a 





Qi 







Q, 





<< 




Q) 


-Qj 





-2L 


2L+ 


L 


Qj 
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H m /a> 


V 







-2L+ 


-H m /m 






Starting with the primary state one can build up a tower of energy eigen- 

states of H ffl by acting the following raising operators (right tower in Fig. 11): 

+ Pi Ho i H co i 

Q] = -J= + iJ -Kt and tf = —--C+-D. (66) 
V2co V 2 2o 2 2 

From the commutation relations [H®, Q\] = ogj and [H a , L7] = 2colJ , Qj or L 1 
acting on |fo) raise its energy eigenvalue by co or 2(0. For example, the states 
(G/)"|^b) with n = 0, 1,2,... have their energy eigenvalues £ = [Aq +n)(d and 
correspond to excitations of the center-of-mass motion in the /-direction, while 
(Ijy^o) have E = [Ao +2n)(0 and correspond to excitations of the breathing 
mode [25]. The primary state \*Fo) is annihilated by the lowering operators Qj and 
L: 

Qi\W o }=0 and L\W o )=0. (67) 

Therefore |fo) corresponds to the bottom of each semi-infinite tower of energy 
eigenstates, which is the ground state with respect to the center-of-mass and breath- 
ing mode excitations. 

The commutation relations of the Hamiltonian and the raising and lowering op- 
erators in the oscillator space are summarized in Table 2. It is clear from the above 
arguments and also from Tables 1 and 2 that the roles of (P,, Kf, D, C, H) in the free 
space is now played by (Qj , Qj, H m , L, L + ) in the oscillator space. Fig. 1 1 illustrates 
the correspondence between the spectrum of scaling dimensions of local operators 
in NRCFT and the energy spectrum of states in a harmonic potential. 

The operator-state correspondence elucidated here allows us to translate the prob- 
lem of finding the energy eigenvalues of the system in a harmonic potential to an- 
other problem of finding the scaling dimensions of primary operators in NRCFT. 
In Sect. 3.3, we use this correspondence to compute the energies of fermions at 
unitarity in a harmonic potential. We note that the similar correspondence between 
quantities in a harmonic potential and in a free space has been discussed in the 
quantum-mechanical language in Refs. [25, 26]. 



Unitary Fermi gas, e expansion, and nonrelativistic conformal field theories 



25 



3.2.3 Unitarity bound of scaling dimensions 

As we mentioned in Sect. 3.2.1, the scaling dimensions of local operators are 
bounded from below. The lower bound is equal to d/2, which can be seen from 
the following intuitive physical argument. According to the operator-state corre- 
spondence, the scaling dimension of a primary operator is the energy eigenvalue of 
particles in a harmonic potential. The latter can be divided into the center of mass 
energy and the energy of the relative motion. The ground state energy of the cen- 
ter of mass motion is (d/2) CO, while the energy of the relative motion has to be 
non-negative [25]. Thus, energy eigenvalues, and hence operator dimensions, are 
bounded from below by d/2. 

More formally, the lower bound can be derived from the requirement of non- 
negative norms of states in our theory [27]. We consider the primary state \*Fo) 
whose mass and energy eigenvalues are given by M -\ and Aq(q: 

M\W )=M ot \To) and H m \W ) = A o (0\W o ). (68) 
We then construct the following state: 

and require that it has a non-negative norm (<f>|<f>) > 0. With the use of the commu- 
tation relations in Table 2 and Eq. (67), the norm of \<P) can be computed as 

(*|*)= Uo-fj (ffo|fo}>0. (70) 

Therefore we find the lower bound on the scaling dimensions of arbitrary primary 
operators: 

&0>\- (71) 

The lower bound d/2 multiplied by co coincides with the ground state energy of a 
single particle in a li-dimensional harmonic potential. 

When the primary state Y¥q) saturates the lower bound Aq =d/2, the vanishing 
norm of {<t>\<P) =0 means that the state itself is identically zero; \<P) =0. Accord- 
ingly, the state created by the corresponding primary operator 0(t, x) obeys the free 
Schrodinger equation: 

O + (f,a;)|0) =0. (72) 

In addition to the trivial one-body operator we will see nontrivial examples of 
primary operators that saturate the lower bound of the scaling dimensions. 

If a theory contains an operator O with its scaling dimension between d/2 and 
(d + 2)/2, then O^O is a relevant deformation: Aq\q < d + 2. Therefore, such a 



2M, 
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theory should contain a fine tuning in the O^O channel. We will see this pattern 
explicitly below. 



3.2.4 Nonuniversality of p-wave resonances 

One consequence of the unitarity bound on operator dimensions is the impossibility 
of achieving universality in p-wave resonances in three spatial dimensions. This 
issue was raised in connection with the a-n scattering [28]. Theoretically, such a 
resonance would be described by the following Lagrangian density: 

£ = v4 (id, + Wo + ci 0, f (id, + ^\ ft + ca^/ft (73) 

+ (yfftv; - ft vtfM- + (ndiVt - 9 ivm) > ( 74 > 

where ci and q are bare couplings chosen to cancel out the divergences in the one- 
loop self-energy of 0,-. (In contrast to the s-wave resonance case, the loop integral 
is cubic divergent and requires two counter terms to regularize, corresponding to 
two simultaneous fine tunings of the scattering length and the effective range. In 
dimensional regularization, c\ = ci = 0.) The field 0, now has a finite propagator 



/ 



-d/2 

dtdxJ**-»"{Tk(t,x)ti(0)) - ( - Po + E--iO+ ) . (75) 



Such a theory might appear healthy but the scaling dimension of 0,-, as one can see 
explicitly by comparing Eq. (75) with the Fourier transform of Eq. (62), is 4^ = 1 
which is below the unitarity bound of d/2 in d = 3. (In a free theory A<f, = d+ 1, but 
the fine tunings "reflect" A§ with respect to [d + 2)/2 so that A^ becomes 1.) Thus, 
p-wave resonances cannot be universal in three spatial dimensions. The proof given 
here is more general than that of Refs. [29] . Other examples considered in Refs. [29] 
can also be analyzed from the light of the unitarity bound. 



3.3 Scaling dimensions of composite operators 

All results derived in Sects. 3.1 and 3.2 can be applied to any NRCFTs. Here we con- 
centrate on the specific system of spin-j fermions in the unitarity limit and study 
various primary operators and their scaling dimensions. The simplest primary op- 
erator is the one-body operator y/ a (x) whose scaling dimension is trivially given 
by 

\o = \- (76) 

This value multiplied by a indeed matches the ground state energy of one fermion 
in a c/-dimensional harmonic potential. 
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The first nontrivial primary operator is the two-body composite operator: (j)(x) = 
c oWiWfi x )' which also appears as an auxiliary field of the Hubbard-Stratonovich 



transformation (8a). The presence of the prefactor cq ~ A guarantees that matrix 
elements of (j)(x) between two states in the Hilbert space are finite. Accordingly its 
scaling dimension becomes 



A$ =A ¥ +A ¥ +(2-d)=2. 



(77) 



This result can be confirmed by computing the two-point correlation function of <p 
and comparing it with the Fourier transform of Eq. (62): 



dtdxe ipo '- ip x (T (f , x)^ (0)) 



d/2 



ri 



-Po- 



2(m t +m^) 



d/2-1 



(78) 

According to the operator-state correspondence, the ground state energy of two 
fermions at unitarity in a harmonic potential is exactly 2co for an arbitrary spatial 
dimension d. This result is consistent with our intuitive pictures of spin-^ and J, 
fermions in the unitarity limit as a single point-like composite boson in d = 4, two 
noninteracting fermions in d = 2, and two identical fermions in d = 1 (see discus- 
sions in Sects. 2.1.1 and 2.3.2). Note that <p(x) in d = 4 saturates the lower bound 
(71) and thus obeys the free Schrodinger equation (72) with mass = m-j- + mj,. 
The same result in d = 3 has been obtained by directly solving the two-body 
Schrodinger equation with a harmonic potential [30], which is consistent with the 
experimental measurement [31]. 



3.3.2 Three-body operators 

We then consider three -body composite operators. The formula to compute their 
scaling dimensions for arbitrary mass ratio m-j- / mj,, angular momentum /, and spatial 
dimension d is derived in Appendix. Here we discuss its physical consequences in 
d = 3. Three-body operators composed of two spin-| and one spin-| fermions with 
orbital angular momentum / = and / = 1 are 

O^ \x)=Z l ^ Vt (x) (79) 

and 

O^ix) =Zf 1 [(mt+mj MVt-'"t(<WVt] ( x )> ( 8 °) 

where i = 1,2,3 and Z/ ~ is the renormalization factor. The mass factors in 
Eq. (80) are necessary so that the operator is primary; [Ki, 0)lI~}' (0)] = 0. 
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Fig. 12 Scaling dimensions of three-body composite operators OiA with angular momentum / = 
(left panel) and / = 1 (right panel) as functions of the mass ratio m-f/wj,. In the right panel, the real 
part (solid curve) and the imaginary part shifted by +2.5 (dashed curve) are plotted. 



(/) 7 

Their scaling dimensions ALh = j + 1 + Yi obtained by solving Eq. (135) are 

plotted in Fig. 12 as functions of the mass ratio rn^Jm^. For / = (left panel), Aj^ 
increases as m-j- //wj, is increased indicating the stronger effective repulsion in the 
s- wave channel. On the other hand, aS for / = 1 (right panel) decreases with 
increasing m^/m^ and eventually, when the mass ratio exceeds the critical value 

m^/m± > 13.607, it becomes complex as A^ 1 ' = | ± ;Tm(yi). 6 In this case, the 
Fourier transform of Eq. (62) implies that the two-point correlation function of 



^! behaves as 



oc sin 



Im(7i)ln 



p 2 — (4f7i| + 2mi)po — 'O 4 
A 1 



(81) 



Now the full scale invariance of the original NRCFT is broken down to a discrete 
scaling symmetry, 



P 



2 WMri) p and Po ^e 2n7l / lm ^p , (82) 



with « being an integer. This is a characteristic of the renormalization-group limit 
cycle and related to the existence of an infinite set of three-body bound states in the 
p-wave channel. Their energy eigenvalues form a geometric spectrum E n+ \/E„ = 
e -2jz/\im(Yi)\ ^ YyjjjQjj j s known as the Efimov effect [32]. Because the system devel- 
ops deep three-body bound states, the corresponding many-body system cannot be 
stable toward collapse. 

We note that, in the range of the mass ratio 8 .6 1 86 < m-j- / < 1 3 .607, the scaling 

dimension of cj^ 1 ' satisfies | < A^J 1 ' < |. Accordingly the following three-body 
interaction term 

£ 3 -body = ci 0$(t,x)0$\(t,x) (83) 



6 This situation illustrates a general feature: onsets of the Efimov effect occur when the ground 
state energy of the corresponding few-body system in a harmonic potential is equal to ff>. 
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becomes renormalizable because now the coupling has the dimension —2 < [c\] < 
0. The Lagrangian density (8a) with £3_b dy added defines a new renormalizable 
theory. In particular, when the coupling c\ is tuned to its nontrivial fixed point (c\ ^ 
0) that describes a three-body resonance in the p-wave channel, the resulting theory 
provides a novel NRCFT. The corresponding system is that of spin-^ fermions with 
both two-body and three-body (ttl) resonances and its many -body physics was 
studied in Ref. [16]. 



3.4 Application of the e expansion 

It would be difficult to determine exact scaling dimensions of composite operators 
with more than three fermions in d = 3. However, it is possible to estimate them with 
the help of the e expansions around d = 4 and d = 2. The formulations developed in 
Sect. 2.2 can be used for the few-body problems as well just by setting ji = 0o = 0. 
Here we concentrate on the equal mass case = mi. 

3.4.1 Scaling dimensions near four spatial dimensions 

The scaling dimensions of composite operators can be determined by studying their 
renormalizations. We start with the simplest three-body operator </>i//| near d = 4, 
which has zero orbital angular momentum / = 0. The leading-order diagram that 
renormalizes (j>y^ is 0(e) and depicted in Fig. 13 (left). Performing the loop in- 
tegration, we find that this diagram is logarithmically divergent at d = 4 and thus 
the renormalized operator differs from the bare one by the renormalization factor; 
(0i//|) ren = •Z^ t ^V / t' where ^ipt = 1 — felnA. From the anomalous dimension 

Y^vt = — InZ^/lnA = f e, we obtain the scaling dimension of the renormalized 
operator up to the next-to-leading order in e: 

=A$+A Vt + y^ =4+-e + 0(e 2 ). (84) 

According to the operator-state correspondence, there is a three-fermion state in a 
harmonic potential with I = and energy equal to = Aq^cq, which continues 
to the first excited state in d = 3. Even within the leading correction in e, the naive 
extrapolation to e — > 1 yields £3°' — > 4.83 (0, which is already close to the true first 
excited state energy 4.66622 (0 in d = 3 [26]. 

The ground state of three fermions in a harmonic potential has / = 1 in d = 3. 
The corresponding primary operator near d = 4 is 2^8^ — (90)y/|, which is 
renormalized by the same diagram in Fig. 13 (left). From its anomalous dimen- 
sion 7= — InZ/lnA = — ^e, we obtain the scaling dimension of the renormalized 
operator up to the next-to-leading order in e: 
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Fig. 13 Leading-order Feynman diagrams to renormalize the composite operators <p tyfr (left) and 
<j><j> (middle) near d = 4 and V'tV-J. (right) near d = 2. These operators are inserted at ®. 



^8v/ t -(8^)v t =^+z\ v , t + l+7=5--e + 0(e 2 ). (85) 

The operator-state correspondence tells us that the three-fermion state in a harmonic 
potential with 1 = 1 has the energy £3 = ^^g^-jg^w CO. The naive extrapolation 

to e — > 1 yields £3 — > 4.17 ft), which is already close to the true ground state energy 
4.27272 CO in d= 3 [26]. 

We now turn to the lowest four-body operator <p<p with zero orbital angular mo- 
mentum / = 0. The first nontrivial correction to its scaling dimension is 0(e 2 ) given 
by the diagram depicted in Fig. 13 (middle). The two-loop integral can be done an- 
alytically and we obtain the renormalization factor; = 1 — 8e 2 In S In A . There- 
fore we find the scaling dimension of the renormalized operator (</>0)ren — ^ths 00 
to be 

AM =2A p +Yw =4 + 8e 2 ln— +0(e 3 ). (86) 

According to the operator-state correspondence, the ground state of four fermions 
in a harmonic potential has the energy = A^^co. Because the 0(e 2 ) correction 
turns out to be large, we shall not directly extrapolate Eq. (86) to e — >• 1 but will 
combine it with an expansion near d = 2. 

The above results can be easily extended to a general number of fermions by 
evaluating the diagrams in Fig. 13 with more boson lines attached. The primary 
operators having N fermion number and orbital angular momentum / are sum- 
marized in Table 3 with their scaling dimensions computed in the e expansion. 
The energy of the corresponding state in a harmonic potential is simply given by 

Ejp = AqCO. The leading-order results [E^ =Nco for even N and E$' = (N+l)co 

and E$ = (N + 2)co for odd N] can be intuitively understood by recalling that 
fermion pairs at unitarity in d = 4 form point-like bosons and they do not inter- 
act with each other or with extra fermions. Therefore the ground state of N = 2n 
fermions consists of n free composite bosons occupying the same lowest energy 
state in a harmonic potential with the energy 2co in d = 4. When N = 2n+ 1, the 
ground state has I = and consists of n composite bosons and one extra fermion 
occupying the same lowest energy state again. In order to create an / = 1 state, one 
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Table 3 Af-body composite operators with angular momentum / near d = 4 and their scaling di- 
mensions in the e expansion. 



N (I) 




A 


2n (I = 0) 
2« + l (/ = 0) 
2n+l (/ = 1) 


<t>" 


N + N(N-2)e 2 In 2? + 0(e 3 ) 
N+ l + ^l e + (e 2 ) 
N + 2+ 2N 1 - 21 e + 0(e 2 ) 



of the n + 1 particles has to be excited to the first excited state, which costs the addi- 
tional energy lo. The leading correction to the energy, which is represented by the 
anomalous dimension in NRCFT, originates from the weak boson-fermion [0(e)] or 
boson-boson [<3(e~)] interaction. Finally we note that in d = 4, we can observe the 

odd-even staggering in the ground state energy; — (£$-1 + £^ +1 ) /2= 1(0 for 
odd N. 



3.4.2 Scaling dimensions near two spatial dimensions 

Similarly, we can determine the scaling dimensions of composite operators near 
d = 2 by studying their renormalizations. Here we consider the three-body opera- 
tors i//|i//j_d f y/j- and y/^y/^dy/i-, which are primary and have the orbital angular mo- 
mentum / = and 1 = 1, respectively. The leading-order diagrams that renormalize 
them are 0(e) and depicted in Fig. 13 (right) with one more fermion line attached. 
Performing the loop integrations, we find that these diagrams are logarithmically 
divergent at d = 2 and thus the renormalized operators differ from the bare ones 
by the renormalization factors; Z^ ¥ g l¥ = 1 + |elnA and Z ¥ ^ ¥ ^g ¥l _ = 1 + |elnA. 

From the anomalous dimensions Jo = — lnZ^/lnA = — |e, we obtain the scaling 
dimensions of the renormalized operators up to the next-to-leading order in e: 

A VtVl d, Vf = + = 5 + 0(e 2 ) (87) 

and 

A Vt¥l a Vl = 3A Va + 1 + 7 rtn a Vt = 4 + 0{l 2 ). (88) 

The operator-state correspondence tells us that the three-fermion states in a har- 
monic potential with / = and / = 1 have the energies £3 = A ¥ . ¥ .^ t ^(0 and 

E^' = A ¥t¥i 9 ¥t (Q, respectively. 

We can develop the same analysis for composite operators with more than three 
fermions by evaluating the diagram in Fig. 13 (right) with more fermion lines at- 
tached. The primary operators having N < 6 fermion number and orbital angular 
momentum / are summarized in Table 4 with their scaling dimensions computed 
in the e expansion. Note that composite operators having the same classical dimen- 
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Table 4 Af-body composite operators with angular momentum / near d = 2 and their scaling di- 
mensions in the e expansion. Known values for the energies of N fermions in a harmonic potential 
in d — 3 are also shown in units of h(0. 



N(l) 


(i) 


A 


E/hm in = 3 


2 {1 = 0) 

3 (/ = 0) 

3 (1=1) 

4 (1 = 0) 

5 (1 = 0) 

5 (1 = 1) 

6 (1 = 0) 


VtYl(dYt) 

(*) 

VtVi(9wr a Vi)9yfr 
VtVi(0Vr 9 V;) 2 


2 

5 + 0(e 2 ) 

4 + 0(e 2 ) 

6- e -+( £ " 2 ) 

9 ll±^05- + o( - 2) 

8-e + 0(e 2 ) 
10-2e + (e 2 ) 


2 [30] 

4.66622 I 26 l 
4.27272 1 26 1 
» 5.028 1 33 1 
« 8.03 1 33 1 
« 7.53 I 33 ! 
^8.48 1 33 1 



(*) = av/ t v/ i (0V / T- s V|)'5Vt + 6Vt 5 iV4.(9'/ t f V4.)^¥t + c^tV|(('3/0V / T)- 9 V / i)'?;VT ~ 
d\ir t Yi(&W dY^idtVt with i<hb,c,d) = (±19^-5^,t16\/3,-6\/35t6v / 3,16\/35). 



sion can mix under the renormalization and thus the primary operator with the well- 
defined scaling dimension may have a complicated form such as for = 5 and / = 0. 
The leading order results for the corresponding energies Ejp = Aq(0 in a harmonic 
potential can be easily understood by recalling that fermions at unitarity become 
noninteracting in d = 2. Therefore the energy eigenvalue of each A^-fermion state is 
just a sum of single particle energies in a harmonic potential in d = 2, and obviously, 
the ground state energy shows the shell structure. The 0(e) correction to the energy, 
which is represented by the anomalous dimension in NRCFT, originates from the 
weak fermion-fermion interaction. We can see in Table 4 the rough agreement of 
the naive extrapolations of Aq to e —¥ 1 with the known values in d = 3. 

3.4.3 Interpolations of e expansions 

We now extract the energy of N fermions in a harmonic potential in d = 3 by in- 
terpolating the two expansions around d = 4 and d = 2 just as has been done for £, 

in Sect. 2.3 and T c /ep in Sect. 2.5. We approximate Ejp /(O by ratios of two poly- 
nomials (Pade approximants) and determine their unknown coefficients so that the 
correct expansions both around d = 4 (Table 3) and d = 2 (Table 4) are reproduced. 

Fig. 14 shows the behaviors of the three-fermion energies EjZL* with orbital an- 
gular momentum / = (left panel) and I = 1 (right panel) as functions of d. The 
middle four curves are the Pade interpolations of the two NLO expansions. Because 
the exact results for arbitrary d can be obtained from Eqs. (135) and (136), we can 
use this case as a benchmark test of our interpolation scheme. We find that the be- 
haviors of the interpolated curves are quite consistent with the exact results even 
within the leading corrections in e and e. In d = 3, these interpolations give 
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Fig. 14 Energies of three fermions in a harmonic potential in the .s-wave channel 1 = (left panel) 
and in the p-wave channel / = 1 (right panel) as functions of the spatial dimension d. The dashed 
(dotted) lines are the extrapolations from the NLO expansions around d = 4 in Eqs. (84) and (85) 
[d = 2 in Eqs. (87) and (88)]. The four solid curves show the Pade interpolations of the two NLO 
expansions. The symbols (x) indicate the exact values for each d obtained from Eqs. (135) and 
(136). 



! 4.71, 4.7, 4.72, 4.72 (89) 



: 4.29, 4.3, 4.32, 4.29, (90) 



E (i=o) 

ft) 

and 

^N=3 
CO 

which span very small intervals E^/co w 4.71 ±0.01 and E^/a 4.30 ±0.02. 
Our interpolated values are reasonably close to the exact results 4.66622 ft) and 
4.27272 ft) in d = 3 [26]. 

Here we comment on the convergence of the e expansions around d = 4 and d = 
2. By performing the expansions up to O(e so ) with the use of the exact formula in 
Eq. (135) and studying their asymptotic behaviors, we can find that the e expansions 
are convergent at least for these three-body problems. Their radii of convergence are 
estimated to be |e| < 0.48 and |e| < 1.0 for the / =0 case and |e| < 1.4 and |e| < 1.0 
for the / = 1 case. 

The same analysis can be done for the energies of more than three fermions in a 
harmonic potential where exact results are not available. The Pade interpolations of 
the two NLO expansions for N = 5 yield 



E (i=o) 

CO 

and 

r 'N=5 
CO 



i 7.71, 7.64, 7.66, 7.82 (91) 



: 7.10, 7.16, 7.19, 7.09 (92) 
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in d = 3, which span relatively small intervals E 5 /co~ 7.73 ±0.09 and E^ /co ps 
7. 14 ± 0.05. On the other hand, the Pade interpolations of the two NLO expansions 
for N = 4, 6 yield 

E ('=o) 

ps 5.55, 4.94, 4.94, 4.90, 6.17 (93) 

co 

and 

E (f=o) 

PS 10.1, 7.92, 7.92, 7.80, 16.4 (94) 

co 

in d = 3. The first and last values corresponding to the Pade approximants where 
all terms are distributed to the numerator or denominator are considerably off from 
the other three values. This would be because of the large NLO corrections near 
d = 4. If such two extreme cases are excluded, the other three values span rather 
small intervals E^/co ps 4.92 ± 0.02 and E^/CO w 7.86 ± 0.06. For comparison, 
the numerical results obtained by using a basis set expansion technique are shown 
in Table 4 [33]. 



4 General coordinate and conformal invariance 

Some nontrivial results can be obtained for the unitary Fermi gas using general 
symmetry arguments. For this end it is convenient to couple the unitary Fermi gas 
to an external gauge field (ji =0,1,2,3) and to an external metric gjj. Both A^ 
and gij can be functions of time and space. Now the action of the unitary Fermi gas 
becomes 

5 = Jdt dx v£ f^- y4 D, y a - l - D, v4 y/ a - ^/ iD i vl D j Wo + w\ + Vt 

(95) 

where D, = d t — iAq and Z), = <9,- — iAj are covariant derivatives. We can see that Ao 
plays the role of the external trapping potential. Recall that when the dimensional 
regularization is used, the term °< 0^0 is absent in the unitarity limit. 

By direct calculations, one can verify that this action is invariant under the fol- 
lowing transformations: 

• Gauge transforms 

f-K*V. (jf^^'^cj) (96a) 
A A - d t a, Aj -t At - diO, (96b) 

• General coordinate transformations 
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x i =x i (t,x i ') (97a) 

yf(t,x) ->• \^(t,x') = y(t,x) (97b) 

Q(t,x)->Q'(t,x') = Q(t,x) (97c) 
<9x' dx- 7 

fty(t,as) ->■ w/fox') = y73-7fty(*,se) (97d) 

A (/,x) ->A,y(*,a;') =A (/,x)+^A i --i f ^y (97e) 

Ai(/,a:) -*A,(t,a/) = -—?Ai(t,x) - —jx J gij (t,x) (97f) 
dx l dx' 

• Conformal transformations 

t->t', t=t(t') (98a) 
( dt\ 3,74 

xif{t,x)^x/(t',x)=(—j \j/(t,x) (98b) 

</.(f,a ; )^0 , (r , , : r)= (Q\+(t,x) (98c) 

Ao(/,x)->A(,(»',aj) = r^Ao(*,x) (98d) 

A i (t,x)->A' i (t',x)=A i (t,x) (98e) 

gy(?,a ; )^4(f',x)= (|^) gi j{t',x). (98f) 



These symmetries allow one to transform the unitary Fermi gas in a free space 
into that in a harmonic potential with an arbitrary time-dependent frequency ©(?). 
This is done by a combination of a conformal transformation t = f(t'), a general 
coordinate transformation 

x i = X- 1 (t)x i ' with A(f) = [/'(f)]- 1 / 2 , (99) 

and a gauge transformation with 

1 A , 

a = -^. (100) 

If one starts with A^ = and gjj = these three transformations leave the gauge 
vector potential A, and the metric gjj unchanged, but generate a scalar potential Ao: 

A = l -co 1 {t)x 2 with co 2 {t) = -^jx 2 . (101) 
The transformed field operator is 
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W'(t,x) =exp ( l -\A X-^{t)W{t),X-\t)x). 




(102) 



This map between the unitary Fermi gas in the free space and that in the harmonic 
potential with the time-dependent frequency was previously found in Ref. [34]. 

In future applications, we only need the infinitesimal forms of the transforma- 
tions. For reference, they are 



for the conformal transformations, where A [O] is the dimension of a field O; A [iff] = 

§, A [0] = 2, A [A ] = 2, A [A,] = 0, and S[ 8ij } = -2. 



4.1 Vanishing bulk viscosities 

One consequence of the general coordinate and conformal invariance is the vanish- 
ing of the bulk viscosity of the unitarity Fermi gas in the normal phase, and the 
vanishing of two (out of three) bulk viscosities in the superfluid phase. These con- 
clusions come from the requirement that hydrodynamic equations describing the 
motion of a fluid in the external gauge field and metric possess the same set of 
symmetries as the microscopic theory. 

4.1.1 Normal phase 

In the normal phase, the hydrodynamic equations are written in term of the local 
mass density p, the local velocity v', and the local entropy per unit mass s. These 
equations are 



<Sy/ 
•SAo 
8A t 



iay-£, k d k y, 

-a-% k d k Ao-A k £ k , 

-djd - ^ k d k Ai-A k d£ k + m gi £ k , 



(103a) 
(103b) 
(103c) 
(103d) 



for the gauge and general coordinate transformations, and 




(104) 




{Ei-F ik v k ), 



d'T 



) 



2R 



(105) 



(106) 



(107) 
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where TIi k is the stress tensor, k is the thermal conductivity, and R is the dissipative 
function. Compared to the usual equations written for the flat metric and in the ab- 
sence of the gauge field, we have replaced the derivatives <3; by the covariant deriva- 
tives Vj and added the force term in the momentum conservation equation (106), 
which comes from the electric force (£, = <9,A,- — <9,Ao) and the magnetic Lorentz 
force (Fx = d{A k — d k Aj). The stress tensor can be written as 

n ik = pvjvjt + pg ik - a' ik , (108) 

where p is the pressure and o' jk is the viscous stress tensor. The information about 
the kinetic coefficients is contained in a' jk and R. 

In the dissipationless limit (a' =R = 0), the hydrodynamic equations are invari- 
ant with respect to the general coordinate transformations, provided that p, s, and v' 
transform as 

8p = (109) 
8s = -$ k d k s, (110) 
Sv'^-^v' + v*^ 1 + «*'•. (Ill) 

Now consider the dissipative terms. To keep the equations consistent with the 
diffeomorphism invariance, one must require that Oy and R transform as a two- 
index tensor and a scalar, respectively: 

dalj = -$ k d k olj - o kj d£ k - a ik d^ k , (112) 
8R = -^ k d k R. (113) 

In a flat space the viscous stress tensor is given by 

a!j = 77 (divj + djvi) + (C - h) 8 iA vk , (1 14 > 

where rj and £ are the shear and bulk viscosities. In the naive extension to the 
curved space where one simply covariantizes the spatial derivatives, a\j is not a 
pure two-index tensor; its variation under the diffeomorphism contains extra terms 
proportional to S, k . The correct extension is 

2 \ r_ k 



°U = ri(V l v J + V i v,+g lJ )+l^--r,) 8ij (V t v*+^ ) • (115) 
Similarly, the dissipative function R in the curved space becomes 



2R = ± \ViVj + V jVi - - 8ij V k v K + gij - - gij °- 



Wivi+ 2gJ + J diTd ' T - (116) 
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We then turn to the conformal invariance. The dissipationless hydrodynamic 
equations are invariant under (104), if the dimensions of different fields are 

A[p]=2A[\jf}=3, A[s]=0, A[v l }=2. (117) 

Now let us consider the dissipation terms. From dimensional analysis, one finds that 
one has to set 

A[t])=A[Q=A[K) = l (118) 

for the hydrodynamic equations to be scale invariant. However, the conformal in- 
variance is not preserved generically. The culprit is gy that transforms as 

8g i j = -Pg i j+'$g iJ , (119) 

which leads to Oy and R not to conform to the pattern of (104), unless the bulk 
viscosity £ vanishes. Thus the requirement of the conformal invariance of the hy- 
drodynamic equations implies £ = 0. 



4.1.2 Superfluid phase 

Similarly, we can repeat the argument for the superfluid case. The hydrodynamics of 
superfluids contains an additional degree of freedom, which is the condensate phase 
6, whose gauge-covariant gradient is the superfluid velocity: 

v^-idiO+At). (120) 
m 

It transforms in the same way as the normal velocity v, = v" under the general coor- 
dinate and conformal transformations. Its consequence is that the relative velocity 
between the superfluid and normal components w' = v' s — v' transforms as a pure 
vector under the diffeomorphism: 

8w' = -^ k d k w' + w k d£'. (121) 

The term in the variation cancels between 8v s and 8 v. 

The diffeomorphism-invariant dissipative function in the curved space is 

2R = 2 { v<Vj + VjVi ~ l 8ijVkX,k + $ ij - \ 8i q 

+ ^ i (p s w')\ 2 + jd i Td'T. (122) 
Under the conformal transformations, R transforms as 
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SR = -PR - l$R + |Ci>3 V,-(pw') + |C2^3 (v.v 1 + JL) . (123) 

The requirement of the conformal invariance of the superfluid hydrodynamics im- 
plies that the j3 terms must have vanishing coefficients, i.e., t,\ = C,% = 0. 

In conclusion, we find that in the unitary limit, the bulk viscosity vanishes in the 
normal phase. In the superfluid phase, two of the three bulk viscosities vanish. 



4.2 Superfluid effective field theory 

At zero temperature, the long-distance dynamics of the unitary Fermi gas is de- 
scribed by an effective field theory, with some effective action S e ff- The effective 
theory should inherit the general coordinate invariance of the microscopic theory. 
This means that the effective action is invariant under the general coordinate trans- 
formations, which in its turn means that 

S£ = -£,%£. (124) 

The low-energy degree of freedom is the phase of the condensate 9. The time rever- 
sal symmetry means that the effective theory is invariant under 

t^f-t and 0->-0. (125) 

Since 9 is a Nambu-Goldstone field, it should always appear with derivatives 
in the effective Lagrangian. Therefore, for power counting purposes, we can set 
9, dj9 ~ 0(p°). The leading-order effective Lagrangian should be 

C = L(9,di9,Ao,Ai,gij)- (126) 

The gauge invariance and the invariance with respect to three-dimensional general 
coordinate transformations (with time-independent E, ') limit the Lagrangian to be a 
function of two variables: 

C = C(D,6, g ij Dj 9Dj9), ( 1 27) 

where 

D,9 = 9-A Q and £>,0 = <?,0 -A,. (128) 

The invariance of the effective theory with respect to general coordinate trans- 
formations requires 

Jj 

C=P(X) with X=D t 9-—Di9D;9. (129) 

2m 

If we now set the metric to be flat and the external field to be zero, we find the most 
general form of the Lagrangian for superfluids: 
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Table 5 Eight classes of scale invariant nonrelativistic systems proposed in Ref. [17]. In all cases 
below, the coupling of the interspecies contact interaction term has the dimension [cq] = — 1 and 
can be tuned to the nontrivial fixed point describing an interspecies resonance. 



# of species and dimensions 


Spatial configurations 


Symmetries + M, D, C, H 


2 species in pure 3D 


x A = x B = (x,y,x) 


Pi, Ki, Jij with i,j = x,y,z 


2 species in 2D-3D mixture 


&a = (x,y) x B = (x,y,z) 


Pi, Kj, Jij with i,j = x,y 


2 species in 1D-3D mixture 


&a = (z) &b = (x,y,z) 


Pz,K : , J xy 


2 species in 2D-2D mixture 


x A = (x,z) x B = (y,z) 


Pz, %z 


2 species in 1D-2D mixture 


x A = (z) x B = (x,y) 


Jxy 


3 species in ID- ID- ID mixture 


x A = 0) x B = (>■) x c = (z) 


None 


3 species in 1D 2 -2D mixture 


x A = x B = (x) x c = (x,y) 


P x , K x 


4 species in pure ID 


x A =x B = x c = x D = (x) 


P.-, K x 



which was previously found by Greiter, Wilczek, and Witten [35] using a differ- 
ent line of arguments. Moreover, by studying the thermodynamics of the effective 
theory one finds that the function P(X) is the same function that determines the 
dependence of the pressure on the chemical potential. 

To the next-to-leading order, the effective Lagrangian contains terms with two 
additional derivatives. The symmetries restrict the number of independent terms in 
the Lagrangian to two. This allows one to relate different physical quantities with 
each other [13]. 



5 Other scale and conformal invariant systems 

We have discussed various theoretical aspects of the unitary Fermi gas. Finally we 
conclude this chapter by introducing other systems exhibiting the nonrelativistic 
scaling and conformal symmetries, to which a part of above results can be applied. 
One such system is a mass-imbalanced Fermi gas with both two-body and three- 
body resonances, which is already mentioned in Sect. 3.3.2. Its many-body physics 
is studied in Ref. [16]. 

The other systems are multi-species Fermi gases in mixed dimensions [17]. In all 
systems listed in Table 5, the coupling of the contact interaction term that involves 
all species has the dimension [cq] = — 1 and thus the theory with such an interaction 
term is renormalizable. In particular, when the coupling co is tuned to its nontrivial 
fixed point (co ^ 0) that describes an interspecies resonance, the resulting system 
becomes scale invariant. We can derive the reduced Schrodinger algebra and the 
operator-state correspondence for such a system [36]. The few -body and many -body 
physics of multi-species Fermi gases in mixed dimensions are studied in Refs. [17, 
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Fig. 15 Feynman diagrams to renormalize three4}ody composite operators. The solid lines are the 
propagators of and y± fields while the dotted lines are the resumed propagators of <j> field. The 
shaded bulbs represent the vertex function Z(p). 



36, 37, 38, 39]. Some of these systems can be in principle realized in ultracold atom 
experiments. Indeed the 2D-3D mixture has been recently realized using 41 K and 
87 Rb and the interspecies scattering resonances were observed [40]. 



6 Appendix: scaling dimensions of three-body operators 

In this Appendix, we derive the formula to compute the scaling dimensions of three- 
body composite operators for arbitrary mass ratio m-j- /m^, angular momentum /, and 
spatial dimension d from a field theory perspective. We first consider a three-body 
operator composed of two spin-^ and one spin-| fermions with zero orbital angular 
momentum / = 0: 

{ ^\x)=Z A l ^( X ), (131) 

where Z^ is a cutoff-dependent renormalization factor. We study the renormaliza- 
tion of the composite operator <j>\ff-[ by evaluating its matrix element (0\<j)\j/^(x)\p 7 —p) . 
Feynman diagrams to renormalize 0m is depicted in Fig. 15. The vertex function 
Z(po,p) in Fig. 15 satisfies the following integral equation: 

z( P0 ,p) = i-ij -^^G^G^-p-m-mk^k) 

= l-j-^G i (-p-k)D(-k)Z(k ,k)\ ko= ^_, (132) 

where we used the analyticity of Z(ko,k) on the lower half plane of ko. G a (p) = 
(po— 2^ + is the fermion propagator and D(p) is the resumed propagator 

2 

of (j) field given in Eq. (78). If we set po = j^, Eq. (132) reduces to the integral 

equation for zip) = Z ( ^_ } p) . 

Because of the scale and rotational invariance of the system, we can assume the 
form of z(p) to be z(p) 06 (x) > where A is a momentum cutoff. Accordingly 
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the renormalization factor becomes Z A °^ A r with y = — d \\\Za /d In A being the 
anomalous dimension of the composite operator 0y/f. In terms of 7, the scaling 

dimension of the renormalized operator is given by 



-r=2 + - + r- 



(133) 



Substituting the expression of z(p) into Eq. (132) and performing the integration 
over |fc| at A — >• 00, we obtain the following equation to determine 7: 



2*V2 



"'i( 2m t+'"i) 



l-d/2 



r(l-f)r(^)sin[(r+l)7r] Jo 



je S in^e Sin[(y+1 ^ (134) 
sinj 



with cos^ = 



cos 9. The integration over 9 can be done analytically in d = 3, 



but otherwise, has to be done numerically. 

Similarly, for general orbital angular momentum /, we can derive the equation 
satisfied by the anomalous dimension 7: 



"|(2mt+ m ;) 
(m|+m^) 2 



l-d/2 



r(l-|)r(^=l)sin[(7/ + / + l)w] Jo 



d9 sin d - 2 9 P, (cos 9 



sin[(yi + l + l)x] 
sinx 



(135) 

where Pi(z) is a Legendre polynomial generalized to d spatial dimensions. 7 The 
scaling dimension of the renormalized operator with orbital angular momen- 
tum / is now given by 



a[J = 4 + A Vt +l + y, = 2 + y + / + 7. 



(136) 



a|1j for / = 0, 1 in d = 3 are plotted as functions of the mass ratio ra| / my in Fig. 12, 



while A^j_ for / = 0, 1 with equal masses = mi are plotted in Fig. 14 as functions 
of the spatial dimension d. 
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